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Multiple Parton Interactions are the tool to obtain information on the correlations between 
partons in the hadron structure. Partons may be correlated in all degrees of freedom and 
all different correlation terms contribute to the cross section. The contributions due to the 
different parton flavors can be isolated, at least to some extent, by selecting properly the 
final state. In the case of high energy proton-proton collisions, the effects of correlations 
in the transverse coordinates and in fractional momenta are, on the contrary, unavoidably 
mixed in the final observables. The standard way to quantify the strength of double parton 
interactions is by the value of the effective cross section and a small value of the effective cross 
section may be originated both by the relatively short transverse distance between the pairs 
of partons undergoing the double interaction and by a large dispersion of the distribution in 
multiplicity of the multi-parton distributions. 

The aim of the present paper is to show how the effects of longitudinal and transverse 
correlations may be disentangled by taking into account the additional information provided 
by double parton interactions in high energy proton-deuteron collisions. 
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1. INTRODUCTION 

Multiple Parton Interactions (MPI) are expected to play an important role at the 
LHC[1]|2]|3]|1]|5][6][7]|8] |9] [TO] p] H2 [□] QI| . While MPI are important to estimate the back- 
ground in various channels of interest for the search of new physics[7j[9|[10j, a further reason of 
interest is that MPI are the tool to obtain information on the multi-parton correlations in the 
hadron structure. The inclusive cross sections of MPI processes depend in fact linearly on the 
multi-parton correlations[15j [16J [17J . 

The experimentally accessible information is represented by the scale factors which character- 
ize the cross section with different numbers of MPI[18j[19j[20j[21J. While all different correlation 
terms contribute to the scale factors, in order to acquire a deeper insight into the non pertur- 
bative properties of the hadron structure, one needs to disentangle the effects of the different 
correlation terms. In the case of high energy proton-proton collisions, the effects of correlations 
of the multi-parton distributions in transverse coordinates and in fractional momenta are all 
mixed in the final observables. The strength of double parton interactions is quantified by the 
value of the effective cross section. The inclusive cross section of double parton interactions in 
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hadron-hadron collisions, ao, is in fact given by 

^D = T - (1) 

2 e ff 

where a a and erg are the inclusive cross sections corresponding to the interactions A and B, m = 1 
when A and -B are identical while m = 2 when the interactions ^4 and B are distinguishable. 
In the simplest case the effective cross section <r e // represents the transverse interaction area 
where double parton interactions take place. A small value of cr e ff may however be originated 
both by the relatively short transverse distance between the pairs of partons undergoing the 
double interaction and by a large dispersion of the distribution in multiplicity of the multi-parton 
distributions. It has been pointed out that the different effects can be separated by studying 
double parton interactions in collisions of hadrons with heavy nuclei[22j[23j. When different 
numbers of target nucleons are involved in the multiple process, longitudinal and transverse 
correlations give in fact significantly different contributions to the cross section. By separating the 
contributions where different numbers of target nucleons are involved in the hard interaction one 
would hence be able to isolate the effects of longitudinal and transverse correlations. Identifying 
unambiguously the final state of MPI may nevertheless represent a serious challenge in the case 
of heavy nuclei, due to the unbalance caused by rescatterings, the energy loss of the recoil jets 
etc. On the other hand dynamics is much simpler and the final state is much cleaner in the case 
of MPI in collisions with light nuclei, which may hence provide a good handle to approach the 
problem. 

The purpose of the present paper is to discuss the simplest case of MPI with a light nucleus, 
namely the case of double parton interactions in proton-deuteron collisions. In our approach 
we will avoid any assumption on the origin of correlations, whereas our goal will rather be to 
enlighten the connections between the different correlation parameters and the physical observ- 
ables. All effects of partonic correlations in fractional momenta and in the transverse coordinates 
will be worked out in detail. The comparison of the double parton scattering cross sections in 
pD and in pp collisions will finally allow a model independent separation of the contributions of 
longitudinal and transverse correlations. 



2. DOUBLE PARTON SCATTERING IN PROTON-DEUTERON INTERACTIONS 
2.1. Both nucleons in the deuteron interact with large momentum exchange 

In Fig. 1 we show the discontinuity of the amplitude A^ 1 representing the contribution to the 
forward proton-deuteron collision amplitude of a double parton scattering diagram where both 
nucleons interact with large momentum exchange. The amplitude is conveniently represented by 
a Feynamn graph, so that the kinematics is correct and the relevant singularities are apparent. 
In the figure the two hard parton-parton interaction amplitudes are denoted T\ and T2, the 
corresponding large pt parton pairs in the final state are pi, P2 and q\, q2, with the positions 
Pi + P2 = P\ Qi + Q2 = Q- The soft vertices (fi p , (j) n represent the (momentum dependent) 
couplings of the interacting partons with the parent proton, neutron and with the corresponding 
low pt final state fragments. When the vertex involves two interacting partons we put a hat, 
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Figure 1: Double parton scattering cross section in proton-deuteron interactions. Both proton and 
neutron interact with large momentun exchange 

(f>. The expressions <&d(D;N) are the relativistic wave function describing the nucleons with 
momenta N and D — N bound in the deuteron; their treatment is discussed in the Appendix. 
The spin variables are not used: we do not consider polarized beams so we can mediate on the spin 
properties, the deuteron is considered a scalar state, the spinorial structure of the components 
disappears in the nonrelativistic limit of the relative motion. 

The flow of momenta in the interaction is shown in the figure. The cross section is given by 
twice the discontinuity of the amplitude of the diagram A^ 1 in Fig.l. The discontinuity is given 
by: 
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,( 2 ) i f JL <h $ 



a2 2 a' 2 

x Ti(l 2 ,a 2 pi,p 2 ) T*(l' 2 ,a 2 ^ pi,p 2 ) T 2 (h,ai -)• 91,92) ^(Z'^ai ->■ 91,92) 

$ D (L>;iV) $* D (D;N') 
X p - iV) 2 - m 2 ][N 2 - m 2 } [{D - N') 2 - m 2 }[N' 2 - m 2 } 
x S(L -h-l 2 - F 3 ) S(L - l[ - l' 2 - F 3 ) 

x S(N - a 2 - F 2 ) 6{N' - a 2 - F 2 ) S(D - N - a ± - Fi) 5{D — N' — a\ - F ± ) 

x 6{h + ai - Q) 5(l[ + a[ - Q) 5{l 2 + a 2 -P) 5(l' 2 + a! 2 - P) 

x Yl d^i d 4 a i d 4 a l i d 4 l i d 4 l l i d 4 F j 5(Fj 2 - M/) d 4 Nd 4 N'd 4 Pd 4 QdMj 2 (2) 

Note that although the expression of DisoA^ 2 ) involves the integration in P± and Q± over all 
the allowed range, we are interested in the regions of the phase space which correspond to a hard 
scattering among partons. The invariant masses of the remnants are Mj, i = 1,2, j = 1,2,3, 
the ten delta-functions express the four-momenta conservations in the vertices, the vectors Fj 
represents the four momenta of the remnants of the nucleons after fragmentation, they will have 
masses Mj certainly larger than the nucleon masses; dQ\ and d£l 2 are the (two-dimensional) 
invariant final state phase space terms of the two hard interactions. 

In the rest frame of the proton-nucleon centre of mass the L + and D_ components grow with 
the cm. energy ^fS while the L_ and D + components decrease as 1/VS. One hence has 

Zi,+, h,+, ai,-, a-2- < y/S 

h ,-, h,-, 01,+, 02,+ < ^ (3) 
Ji,t, h,ti ai,t, «2,t ^ j=> 

being R of the order of the nucleon size and independent of S, the same holds for the I' and a' 
variables. 

By making use of the delta-functions we may write 



h = (L - F 3 )/2 + A l[ = (L - F 3 )/2 + A' 

l 2 = (L- F 3 )/2 - A l' 2 = (L- F 3 )/2 - X' 

ai = {D-Fi- F 2 )/2 + a a' 1 = {D-F l - F 2 )/2 + a' 

a 2 = \d — F\ — F 2 )/2 -a a 2 = (D — F± — F 2 )/2 - a' 

N = (D-F 1 + F 2 )/2 -a N' = (D-F± + F 2 )/2 - a (4) 

which, besides the overall four-vectorial conservation L + D = Fj + P + Q, leaves us with the 
integrations over d 4 X, d 4 X', d 4 a, d 4 a', dFj, d 4 P, d 4 Q and with the remaining delta-functions [1] 



[1] apparently there is a missing 8— function: had we chosen {L, D} 7^ {!/, D'} we should have found a 8(L + D — 
L' — D'), in the actual case we factor out the <5(0). 
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4<S(A + a + (P - Q)/2)5(X' + a' + (P - Q)/2)S(L + D - ^ F, - P - Q). 

j 

To perform the integration over the longitudinal degrees of freedom, we notice that in the 
upper part of the diagram the minus— components are small while the plus— components are 
small in the lower part of the diagram. Hence the integration on X-,X'_, whose range of 
variation is of 0(1/VS), affect only variables with an equally small range of variation of the 
minus— components, i.e. in the upper part of the diagram, while the minus— components in 
the lower part of the diagram, whose range of variation is of O(VS), remain unaffected. The 
opposite property holds for all variables of type a+, a' + . 

At fixed values of Fj , P, Q the integration in A_ , X'_ is hence performed by disregarding the 
delta-functions above, since in the delta-functions all minus— variables, but A_, are large. The 
same attitude is taken when integrating in a+,a' + . The integrations in A+, A+, ct-, ot'_ are 
performed, on the contrary, by using the delta-functions, where all the small terms are neglected. 
In so doing one obtains A+ = A' + = (Q — P) + /2 and a_ = a!_ = (Q — P)_/2, which in turn 
implies A_ = N'_, h+ = Q+, h+ = P+, «i- = Q-, a,2- = P-- 

We proceed by defining the two-parton amplitude in the nucleon as 

, = 1 [fa d\- 1 [ fa d\'_ 

W2 V2J 2iri V2 l'^ 2-Ki ' 

The integration over a + involves the vertices (f) p ,(f) n ,^D and the denominator 



[(£> _ TV) 2 - m 2 } ■ [N 2 - m 2 ] -a 2 -a 2 

In the variable a + there are then four simple poles, two with positive and two with negative 
imaginary part. We can consider the two poles originating from the request a 2 . = 0, or A^ 2 = m 2 . 
The first choice forces iV 2 to be at least as large as M 2 , so the nucleon is strongly off mass shell 
as compared with the binding energy of the deuteron. The pole N 2 = m 2 implies a space-like 
02, which can however be near a| = 0. The integrations in a + ,a' + are hence approximated by 
keeping only the contribution of the nucleon poles. In this way we can define, using the relation 
(Eq. 3) between and a, the covariant amplitude for finding a nucleon in the deuteron 



If ® D da + _ 1 1 ® D 

^p2j [{D - Nf - m 2 ] ■ [N 2 - m 2 ] 2wi ~ [(D - N) 2 - m 2 ] 

and in strict analogy we have the twin relation: 



D 



da'. 



1 1 



D 



y/2J [(D - A') 2 - m 2 ] ■ [A' 2 - m 2 } 2ni ^2 N'_ [(D - A') 2 - m 2 ] 



N' 2 -- 



A' 



(6) 



In this way the denominators a 2 are left as factors, in analogy with a previous definition we 



set 



for the one-parton amplitude one must here remember that the variables are, as functions 
of a, subjected the condition coming from N 2 = m 2 . The dependence of the function ipn and 
ip2 on the invariant mass of the residual hadron fragments is understood. 

The integration of the longitudinal momenta of the remnants, Fj, is performed by using the 
equivalence d 4 F5 + (F 2 - M 2 ) = \dF±/F± x d 2 F t . One may notice that F%_ is very small and 
-^1+ ^2+ are very small too. The conditions expressed by S(L + D — ^2j ,Pj~P~Q) can hence be 
used to define -F3+, disregarding the effect of (Fi+F 2 ) + , and to define (F1+F2) — , disregarding the 
effect of F 3 _. In so doing from dF 2 _/F 2 _, dF l _/F 1 _ one gets dK-/[(L + D - P- Q) 2 _/A - K 2 ], 
where K_ = (F 2 - Fi)_/2 and it is related to iV_ by the relation: K- = (N + Q - P - D/2)_. 
At fixed P_, Q_ hence dK_ = dN_. 

We can now proceed with the integration on the transverse variables, in the frame where 
both Lf, and Df, are equal to zero. We take now the two-dimensional Fourier transforms (it is 
understood that all the variables are two-dimensional vectors). 



and analogously for the complex conjugated functions, with the variables b^, b' 2 , /3[, f3' 2 , B' . 

The integration over the transverse-momentum variables gives the diagonal property b\ = b^ 
and so on. Moreover, as shown in Fig. 2, one obtains the geometrical relation: 61 — 62 = B — fi\—^ 2 - 

It is useful to introduce now the fractional plus or minus momenta in the following way: 



We introduce also the fractional momentum of the nucleon inside the deuteron Z = 2N-/D- 
so that the fractional momenta of the partons with respect to their parent nucleons are z\ = 



The hard scattering amplitudes are treated as elastic scatterig between massless partons. 
Hadronization is not included. When neglecting the dependence on kinematical factors which 
do not grow with S one obtains: 




(7) 



xi = h+/L+, x 2 = h+/L+, x[ = 2ai_/D_, x' 2 



2a 2 /D 



x' 1 /(2 - Z), z 2 = x' 2 /Z. 



1 



T 1 \ 2 8 A {P -pi-p 2 )5{p 2 1 )5{p 2 2 )d 4 pid A p 2 = 2Sx 1 x' 1 da(xix' 1 ,p 1 j) 



(27T) 2 



1 



T 2 \ 2 5\Q - qi - q 2 )5{q 2 )5{q 2 )d\ 1 d i 



q 2 = 2Sx 2 x' 2 da(x 2 x 2 ,q ltt ) 



(8) 



(2n) 2 



In term of the variables Xi, Z one has: 
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Figure 2: Transverse coordinates in the double parton collision with two target nucleons involved 



dK- = dN- = (D-/2)dZ 

F 3+ = L + (l - xi - x 2 ) 

F x _ = (D_/2)(2-Z)(l-Z!) 

F 2 _ = (D_/2)Z(1 - z 2 ) • (9) 

The one-body and two-body parton densities are defined by the following integrals on the 
invariant mass of the residual hadron fragments: 

T(z,f3) = [ |^ M ( z ,/3)| 2 -^dM 2 



2(2tt)3 J i 

1 f \.l._r„ .. i. lm2 2P +Q+ , !U - 



r(x 1 ,x 2 ,b 1 ,b 2 ) = — — g / \i> M (xi,x 2 ,b 1 ,b 2 )\ 2 - — ^^dM 2 (io) 

4(271-)° J 1 — xi — x 2 

The dependence on pi tQl,t can be transformed into an angular dependence on Oi, and, finally, 
dP±dQ± = S 2 dxidx 2 dx' 1 dx' 2 - The final expression of the double parton scattering cross section 
where both nucleons interact cr^^J is hence given by 

a doubie\ 2 =7^)3 j^(x u x 2 ,b 1 ,b 2 )-^--^-T(x' 1 /Z,l3 1 )r(x , 2 /(2-Z),l3 2 ) 

\^> D (Z,B)\ 2 dBdb l db 2 dPx d(3 2 5{B + b 1 -b 2 - A + /3 2 ) 

[Z(2 - Z)]~ 1 darida:2«ic' 1 d4dZdnidJ22 (H) 

where, as discussed in the Appendix, ^o{Z,B) is expressed in terms of light cone variables 
through the non relativistic deuteron wave function (fp(p 2 ) as 
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and (as already mentioned in the Appendix) the expression has to be understood as being 
symmetrized for the exchange Z <^ (2 — Z). Notice that, since 

J\^ D (Z,B)\ 2 dZd 2 B = J \y D (Z, Pt )\ 2 dZd 2 pt, \^ D (Z, Pt )\ 2 = (2n) 3 E p \<p P (p)\ 2 (13) 

and, as readily verified by using the explicit expressions of E p and p z in the Appendix, dZ/Z = 
dp z /Ep, one has 

^3 j \m D (Z,p t )\ 2d ^d 2 p t = f \vp(p)\ 2 d 3 p = 1 (14) 

Given the symmetry of ^f£>(Z,pt) for the exchange of Z with 2 — Z one moreover can easily 
prove the relations 



(2vr) 3 

which show that all densities in Eq.(ll) are properly normalized. 

2.2. Only one nucleon in the deuteron interacts with large momentum exchange 

In Fig. 2 we show the discontinuity of the amplitude representing the contribution to the 
forward proton-deuteron collision amplitude of a double parton scattering diagram where one 
nucleon interacts with large momentum exchange and the other is a spectator. 
The expression of the discontinuity is 



DisaAW - — f 

A ~ (2tt) 21 J 



7 2 7 2 7/ 27/ 2 ~2~2 /2 /2 
<1 «2 1 1 l 2 a l a 2 0,\ «2 



x Ti(l 2 ,a 2 ->Pi,P2) T*(l' 2 ,a 2 -> pi,P2) T 2 (h,a\ -)• 91,92) T 2 *(Z' 1 ,a / 1 ->■ 91,92) 



[(£> - N) 2 - m 2 ] [{D - N') 2 - m 2 ] 
x 5(L -h-l 2 - F 3 ) S(L - l[ - 1' 2 - F 3 ) 

x 6(D - N - ai - a 2 - F ± ) 5{D — N' — a[ - a 2 - Fi)<5(JV - F 2 )6{N' - F 2 ) 

x 5(h +ai-Q) 5(l[ + a[ - Q) 5{l 2 + a 2 -P) 6{l' 2 + a! 2 - P) 

x Y[ dn { d 4 a i d 4 a' i d 4 l i d 4 l' i d 4 F j 5(F j 2 - Mj 2 ) d 4 Nd 4 N'd 4 Pd 4 Q (16) 
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l 2 '=P + +k/2-8 



a,'=P -k/2+8'+P 



l,'=Q + +k/2+8' 




N / a| '=Q.-k/2-8'+Q t 



Figure 3: Double parton scattering cross section in proton-deuteron interactions. The neutron is a 
spectator in the large momentum exchange process 

Also in this case, beyond yielding the overall four- vectorial conservation L-\-D = Fj+P+Q, 
the delta-functions allow us to write 



h 
h 
a i 

«2 



(L-F 3 )/2 + A 
(L-F 3 )/2-X 
(D-Fi- F 2 )/2 + a 
(D-Fx- F 2 )/2 



a 



[L — F 3 )/2 + A' 

(L-F 3 )/2-A' 
= (£) - Ft - F 2 )/2 + a' 
= (D-F 1 - F 2 )/2 



a 



117) 



and obviously N = F 2 = N' . Multiple integrations are performed analogously to the 
previous paragraph. One needs to introduce the two-parton wave-function also for the bound 
nucleon of the deuteron (the one defined carrying four-momentum S in Fig. 3). So we define the 
two-parton amplitudes of the bound nucleon as 
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V'2 



1 f <t> P da+ _J_[_fp da '~ 

V2j a 2 a 2 2iti ^ ^2 J a[ 



+ 



a' 2 a' 2 2vri 



The spectator nucleoli is set on mass shell. Notice that, in spite of that, one may still claim 
that final state interactions of the spectator are taken into account. The statement is supported 
by unitarity: If the nucleon is produced on mass shell and undergoes a final state interaction 
with the remnants of the other nucleon, final state interaction does not modify the inclusive 
cross section, since the spectator is not observed. If the nucleon is produced off mass shell, its 
virtuality is anyhow rather small and it may not be unreasonable to extend the unitarity relation 
SS^ = 1 to the actual kinematical domain. Thus unitarity allows us to replace the whole final 
state with a cut nucleon line, i.e. with a nucleon on mass shell. 

The different conservation relations, with respect to the case discussed in the previous para- 
graph, give different relations in the transverse coordinates and one obtains that the integration 
on the B coordiante is decoupled from the others which, on the contrary, are linked by the 
relation b\ — b 2 = ft — ft ■ 

By following the steps and introducing the functions and the variables of the previous para- 
graph and keeping into account the presence of two target nucleons, one obtains the expression 



"Swell =z^z j \MP+,Q+MM)mA 2 zp + P-^Q + Q- 



da 



dQ 2 

x\ip 2 {P-,Q-;l3i,/32)M 2 \' 2 \^D(Z,B)\ 2 dBdb 1 db 2 dfa d/3 2 8(b 1 - b 2 - ft + ft) 

dP + dP-dQ + dQ-dN- m tn - , « 1 . . 

rffli dn 2 dM?dM$ (18) 



iV_(iV_ -Q_ -P_)(L+-P + -g+) 1 "8(2vr) 12 

The double parton distribution of the bound nucleon is given by 

1 If 2P-0- 

-r(xi/z,4/z ; ft,ft) = J I ^m 2 (P- , Q- ; h , b 2 ) | 2 - _ dM 2 (19) 

and the cross section is readily expressed as: 

° P d ouUe\l = (^/ r ^'^,6l,6 2 )^-J|r(x' 1 /Z,x' 2 /Z,ft,ft) 

\V D (Z, B)\ 2 [Z{2 - Z)\~ l dB dbi db 2 dft dft S(h - b 2 - ft + ft) 
dx\dx 2 dx\dx' 2 dZ dVL x dQ 2 (20) 

Notice that, as an effect of the nucleon motion, \^d(Z, B)\ 2 is coupled to the interactions by 
the integration on the fractional momentum Z, while the integration on the transverse variable 
B is decoupled from the other transverse variables. 



3. SOME GENERAL REMARKS ON CORRELATIONS 



To the purpose of gaining some understanding of the effect of the correlations in the hadron 
structure on the double parton scattering cross section, in the present paragraph we discuss some 
basic properties of the correlation functions. 
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The expression represents the one-body density of a set of identical particles, as a function 
of the variable £. The density is normalized to the average multiplicity (n) 

Jp(Z)d£ = (n). (21) 

The two-body density p(£i,£2) is analogously normalized to the second moment of the mul- 
tiplicity distribution and, in general the n-body density is normalized to the nth moment of the 
multiplicity distribution 

J pfa... d£ n = {n-(n-l)... 2-1) (22) 

Correlations are introduced in the double parton distribution as deviations from the Poissonian. 
One may hence express 



(n(n — 1)) 



Pfo'h) = - p(aM6) + r?(6,6) (23) 



where the term 77(^1, £2) is such that 

J r/(6,6Rid6 = (24) 

In the uncorrelated case p(£i,£ 2 ) = p{€i)p{&) an d one has (n(n — 1)) = (n) 2 , 77(^1,^2) = 0. The 
presence of correlations thus imply (n(n — 1)) 7^ (n) 2 and/or 7/(^1,^2) ^ 0. In such a case, the 
distribution in multiplicity is not a Poissonian or the dependence on £1, £2 is not factorizable as 
a product p(£i)p(£ 2 )j or both. 

For the double parton distribution one may hence write 



r(xi, x 2 ; 61, 6 2 ) = K(x 1 ,x 2 )F(x 1 ,b 1 )T(x 2 , b 2 ) + C(xi, x 2 ; h, b 2 ) (25) 

where 



if(xi,ar 2 ) - 



G(si)G(s 2 ) 
G(x) = y T(x,6)d6 

G(xi,x 2 ) = J T(xi,x 2 ;bi,b 2 )dbidb 2 (26) 

and 

J C(x 1 ,x 2 ;b 1 ,b 2 )db 1 db 2 = (27) 
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All information on correlations is therefore contained in the two terms K(xi,x 2 ) and 
C(xi, X2\ bi, 62). When K = 1 partons are uncorrelated after integrating in the transverse 
variables, while all information on the correlations on the transverse variables is contained in 
C(x 1 ,x 2 ;bi,b 2 ). 

The expression of the double parton scattering cross section in proton-proton collisions may 
hence be written as a sum of four terms 

m r 2 

a doubie = -J T(xi,X2]b l ,b2)a(xi,x[)a(x 2 ,x' 2 )T(x[,x' 2 -,b l - B,b 2 - B)Y\dx i dx' i db i ddB 
^ i=i 

Tfl f ( 

= — ^K(x 1 ,x 2 )T(x 1 ,b 1 )r(x2,b 2 )a(x 1 ,x' 1 ) 

xd(x 2 , x' 2 )K{x\, x / 2 )r(x / 1 ,6i - B)Y{x' 2 M - B) 
+K(x 1 ,x 2 )T(x 1 ,b 1 )r(x 2 , b 2 )a(xi,x' 1 ) 

xd(x2, x' 2 )C(x' l , x' 2 ; bi - B,b 2 - B) 
+C(xi, x 2 ; bi,b 2 )a(xi,x' 1 ) 

xa(x2,X2)K(x' l ,X2)r(x' 1 ,b l - B)T(x 2 , b 2 - B) 
+C(xi, x 2 ; b 1 ,b 2 )a(xi,x' 1 ) 

2 

xa(x 2 ,X2)C(x' 1 ,x , 2;b 1 -B,b 2 - £)} JJ dxidx^dB (28) 

where m = 1 in the case of two identical hard interactions and m = 2 in the case of two 
different hard interactions. The hard parton-parton cross sections a are integrated with a cutoff 
in momentum transfer. The four contributions to the cross section are represented in Fig. 4. 

As evident in Fig. 4 the unknown quantities related to the different source of correlations, in 
fractional momenta and in the transverse coordinates, cannot be disentangled by measuring the 
double parton scattering cross section only in proton-proton interactions. Additional independent 
information may be however obtained by measuring double parton collisions in proton deuteron 
interactions. The different contributions to the cross section are shown in Fig. 5 

As apparent looking at the two figures 4 and 5, by comparing the double parton scattering 
cross sections in pp and in pD interactions one has the possibility to decouple the effects of longitu- 
dinal and transverse parton correlations, represented by the terms K(x\,X2) and C(xi, X2\ ^1,^2) 
in the double parton distributions. 

4. A SIMPLE ESTIMATE OF THE DOUBLE PARTON CROSS SECTION 
4.1. Effects of the correlations in the transverse coordinates 

Present (still limited) evidence of double parton collisions is consistent with an effective 
cross section independent on fractional momenta. It wouldn't hence be inconsistent to assume 
K(xi,X2) = K, constant. A simplest expression of the double parton distribution where all 
properties previously discussed are implemented is the following factorized expressions 
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Figure 4: Different contributions to the double parton scattering cross section in proton-proton collisions 



T(x;b) = G(x)- 



-b 2 /R 2 



irR 2 

e -(62+&!-Abi-b 2 )/(K 2 [l-A 2 /4]) 
(1 - A 2 /4)(vri? 2 ) 2 



r(xi,x 2 ; 61,62) = KG{x x )G{x 2 ) ^ ^ /AW _ n9 , 2 (29) 



which satisfies 



V{x ll X2-MM)db 2 = KG{x l )G{x 2 ) e -^- = KG{x 2 )T(x 1 ,b 1 ) (30) 

The two parameters K and A represent two different sources of correlations: K 7^ 1 implies a 
non Poissonian multiplicity distribution of partons in the relevant range of fractional momenta, 
while A / implies a non trivial correlations of partons in the relative transverse coordinates. 
The resulting expression of the double parton scattering cross section, in the case of identical 
interactions in pp collisions, is 



pp 

double 



1 r , 2 

- / T(xi,x 2 ; 61, b 2 )a(xi,x , 1 )a(x 2 ,x 2 )r(x , 1 ,X2; h - B, b 2 - B) dxidx\dbidB 

i=l 

1 K 2 a% 

2 2vri? 2 (2-A) ^ 
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Figure 5: Different contributions to the double parton scattering cross section in proton-deuteron collisions 



where o g is the single scattering cross section of the QCD parton model. The effective cross 
section is hence related in a simple way to the correlation parameters K and A 



_ 27^(2 -A) 

The small observed value of the effective cross section may hence imply non trivial correlations of 
partons in the relative transverse coordinates (A > 0) or a non Poissonian multiplicity distribution 
of partons at small fractional momenta {K > 1) or both. This simplest model also shows that 
the two sources of correlation cannot be disentangled by measuring double parton collisions in pp 
interactions only. An independent constraint is however provided by measuring double parton 
collisions in proton-nucleus interactions. The double parton cross section with a spectator nucleon 
in fact is given by 

ill = j^jn^^MM^^-n^/z^yz,^,^) 

\*d(Z, B)\ 2 [Z{2 - Z)]- l dB dh db 2 dtPx dp 2 5(h -b 2 -Pi + /3 2 ) 
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dx\dx 2 dx' x dx 2 dZ dQ.\ dfl 2 



2K 2 



2 _ R2{2 A) " J as(x 1 x' 1 /Z)as(x2x' 2 /Z)\^ D (Z,/3)\ 2 [(27TfZ(2 - Z))~ l dp dZ 

= 2 / aZ uUe (x 1 ,x[/Z,x 2 ,xyz)\^ D (Z^)\ 2 [(2nfZ(2-Z)]- 1 d^dZ 

- 2<ou b ie(*ux[,x 2 ,x 2 ) (33) 

while the double parton cross section where both nucleons take part to the double hard interaction 
is given by 



double 1 2 



T(x 1 ,x 2 ,b 1 ,b 2 )^-^-T(x' 1 /Z,f3 1 )T(x' 2 /(2-Z),(3 2 ) 

\y D {Z,B)\ 2 dBdh db 2 dfc dp 2 5{B + h - b 2 - ft + fo) 
[(2vr) 3 Z(2 - Z)] _1 dxxdx 2 dx' x dx' 2 dZ d£l x dSl 2 

K J as(x 1 x' 1 /Z)as(x 2 x' 2 /(2-Z)) 

xg(p, R 2 (4 - X))\if D (Z, 0)\ 2 [(2vr) 3 Z(2 - Z)\ ^dpdZ 

K j as{x 1 x' 1 /Z)as{x 2 x' 2 /(2-Z))\^ D (Z,0)\ 2 [(27rfZ(2-Z) 

K °eff J aZ uble (x l ,x' l /Z,x 2 ,x' 2 /(2-Z))\^ D (Z,0)\ 2 [(^fZ(2-Z)}- 1 dZ 



dZ 



K(T eff 

ttRI 



(34) 



where the function g(/3,i? 2 (4 — A)) is a Gaussian normalized to one, with argument j3 and 
radius [it! 2 (4 — A)] 1 / 2 . The motion of the nucleons inside the deuteron is slow, their frac- 
tional longitudinal momentum Z is approximately 1, so the nucleons distribution D(Z,B) = 
\^ d (Z,B)\ 2 [(2tt) 3 Z(2 - Z)y X has been treated as a function 5(1 - Z). 

The term o^ouMeL re P resen t s the contribution to the pD cross section with additional in- 
formation on the correlations. By comparing cr^ouMeL with the double parton scattering cross 
section in pp interactions one may have an indication on the value of K which, measured at dif- 
ferent values of x±, x 2 , allows to access the information on longitudinal correlations in a model 
independent way. As an order of magnitude estimate the contribution of the term c^ Me | 2 to 
the proton deuteron double parton scattering cross section is of order a e ff / '(2irR 2 D ) ~ 5%. 

4.2. Smearing with the Deuteron wave function along the longitudinal direction 

In the previous subsection the nucleons distribution D(Z,B) has been treated as a function 
(5(1 — Z) when integrating on Z. One may estimate the effect of the finite width of the deuteron 
by expanding the integrand around 1 as follows 



-> J F(x/Z)D(Z)dZ, F( X/ Z ) = F( XH 



dF 
dZ 



d 2 F 



(Z - If 



Z=l 



(35) 
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Given the symmetry of D{Z) around 1, the leading correction to the (5-function contribution 
is due to the second derivative term, which is hence multiplied by \ f(Z — l) 2 D(Z)dZ. By 
expanding the expression of p 2 in the appendix near Z = 1 one obtains 

p 2 = p 2 + ^(Z-l) 2 + ... (36) 

In the Deuteron wave function (Z — l) 2 is hence multiplied by i? 2 M|,/4, with R a length of 
the order of the radius of the Deuteron. One may hence estimate 



1 -J(Z-lfD(Z)dZ^ J ^^0.5% 



(37) 



The effect is hence roughly one order of magnitude smaller as compared to the correction due to 
the transverse structure of the Deuteron. 



5. CONCLUSIONS 

The inclusive cross sections of multiple parton interactions depend linearly on the multiparton 
correlations of the hadron structure. Partons may however be correlated in all their different 
degrees of freedom, while all different correlation terms contribute to the cross section. 

Whereas correlations in the flavor indices may be studied by selecting properly the final states, 
in this paper we have focused on the problem of disentangling the contributions of the correlations 
in the transverse parton coordinates from the correlations in the distribution in multiplicity of the 
multi-parton distributions. In the case of pp interactions the two contributions are unavoidably 
mixed. The double parton scattering cross section in pD collisions can however provide additional 
information, which may be used to decouple the two terms. 

In the present paper we have worked out in detail the two different contributions to the dou- 
ble parton scattering cross section in pD interactions, corresponding to the different possibilities 
of having 1) just a single nucleon or 2) both nucleons taking part to the hard process. In the 
first instance the integration on the relative transverse coordinate of the two nucleons is decou- 
pled from the other transverse variables and the scale factor in the corresponding cross section, 
a double] v * s ( a P ar t from minor correction terms) the same which one finds in nucleon-nucleon 
collisions, namely cr e ff- As a consequence c^^g^ is twice as big as the double parton scattering 
cross section in pp collisions. In the second instance the scale factor, in the corresponding cross 
section c^ouMeL' is mainly provided by the transverse size of the deuteron. As a consequence the 
value of o'douMeL * s ver y sensitive to the actual value of the correlation term K(x\,X2), which 
measures how much the multi-parton distribution in multiplicity is different from a Poissonian, 
at given values of fractional momenta. 

The first term, o^ouMeli' is the term usually taken into account when considering large p t 
processes in hadron-nucleus collisions. The second term, c^ou&ZeL' represents a positive (namely 
a antishadowing) correction. A rough quantitative estimate indicates that cr^ uble \ 2 may be of 
the order of 5% of Odou&Zelr Apart from the terms representing the double parton correlation, 
all inputs are well known quantities. One may hence conclude that, by measuring double parton 
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collisions in pD interactions, even a reasonably limited statistics should allow obtaining a reliable 
indication on the size of K(x\,X2) which, in turn, will permit a model independent estimate of 
the typical transverse distances between different parton pairs of the structure of the hadron. 

The present discussions shows that pD and, more in general, the interactions of hadrons with 
light nuclei represent the tool to measure the multi-parton correlations of the hadron structure, 
allowing not only to measure the scale factors, which quantify the rate of the different multiple 
parton interactions in proton-proton collisions, but also to obtain, for the first time, the additional 
information needed to disentangle the different correlation terms of the multiparton structure of 
the hadron. Although interactions of deuterons with nuclei and between polarized proton beams 
have been studied at RHIC, no plans are unfortunately foreseen at RHIC to study the interactions 
of protons with light nuclei. In spite of the limited rate of double parton collisions expected 
at the RHIC energy, double parton collisions have been measured by the AFS collaboration 
in proton-proton collisions at the CERN ISR[18j, where the cm. energy was sizably smaller. 
Measuremets of double parton collisions in pD interactions should hence be feasible at RHIC, 
while the inclusion of the study of pD interactions in the RHIC Physics Program would give the 
possibility to obtain unprecedented information on the two-body correlations between partons 
of the hadron structure and hence on the three dimensional parton structure of the hadron. Of 
course the study of double parton collisions in pD interactions would be greatly facilitated at 
the LHC, where the much larger parton luminosity would allow to perform a detailed study of 
various multi-parton scattering channels in a much broader range of fractional momenta, which 
might represents a good motivation to support the option of running, at some stage, light nuclear 
beams at the LHC. 



Appendix A: The deuteron wave function 



One may consider the coupling of the photon to the deuteron as shown in fig 6. 
In the simplest case the vertex is a constant and the nucleons can be treated as scalars. 
The expression of the diagram is 



2Dn - 2p^ + ^ 



D 



d 4 p 



(D — p) 2 — m 2 + ie p 2 — m 2 + ie (D — p + q) 2 — m 2 + ie (27r) 4 



2D^ + Qfl (Al) 



By taking the '0' component in the deuteron rest frame and going to the limit q — > one must 
obtain: 



M D - po dp Q d?p 



Po 



.(M D -p ) 2 -E 2 + ie 
where E p = \/(p 2 + m 2 ). An equivalent expression is 



E 2 + ie M D 



(A2) 







2M D dM D 



(M D - po) 2 - E 2 + ielp 2 -E 2 + ie (2tt) 



dpod 3 p 



(A3) 
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Figure 6: Photon-deuteron vertex 



One may integrate by taking the residues at po = E p — ie and at po = E p + M/j — ie 



-2tt 
2M d 

-1 
2M D 



d 
dM D 

d 
dM D 
-d 



1 



E 2 

v 



dM D 



2E P (Mb - E p ) 2 - 
1 2M 2 D § 2 D 
2E~M D - (2E p M D ) 2 \ (2vr) 3 
d 3 y> 



J 

2E P (M D + E p ) 2 
d 3 p 



E 2 



d 3 p 
j2nT 



2& D 



Ml 



±E 2 



2E p (2n) 3 

one hence obtains the normalization condition 



2& D 



l(M 2 D -AE 2 ) 2 \2E p (27T) 



d 3 p 



(A4) 



Notice that 



2<S> 2 D d 3 p 
(M 2 d -4E 2 ) 2 2E p (2ttY 



(A5) 



D 



{D-pf 



4£ 2 



-4(ra5 + p 2 ) 



(A6) 



where the last equation holds in the non relativistic limit and B is the deuteron binding energy 
(Mj3 = 2m — B). For constant the last expression above coincides with the asymptotic limit 
at large distances of the deuteron wave function: 
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More in general one may obtain a relation between the vertex function <£> d and the non relativistic 
solution of the Schrodinger equation by considering the solution of the Bethe-Salpeter equation 
in the scalar case, in ladder approximation j2_4j [2 5j 



x{p,p) 



f +P 



m 2 + ie 



P 



m 2 + ie 



;x(p,p + q) 



(A8) 



P/2+p 




P/2-p 



Figure 7: Bethe-Salpeter equation 



where P is the deuteron momentum and p is the nucleon momentum in the deuteron rest frame 
(see Fig. 7). The binding is given by the exchange of a particle of mass fx coupled to the nucleons 
with the constant g. The vertex function <&_d is hence expressed in terms of the solution of the 
Bethe-Salpeter equation x as: 

*d(p) = [ -^- 2 x(P,P + q)^i (A9) 
J q 2 - fj, 2 (2tt) 4 

In the deuteron rest frame one assumes instantaneous interaction, in such a way that 



ig 2 ~ig 2 , f ig 2 /rt N d 4 q f -ig 2 , . d 3 q 

q A - fi 2 + ^ J q - M (2vr) 4 J + \i* (2iry 

where 



^p(p + q) = I x(P,p + q)^ 



(All) 
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In the deuteron rest frame one may hence integrate the Bethe-Salpeter equation in qo as follows 



/ 



x(p,p) 



dqp 
2ir 



( P V 



-ig- 



m + it 

d 3 q 



P 



p I — m + it 



dqp 
2ir 



2 +/i 2 



<pp(p + q) 



(2tt) ; 



(A12) 



which gives 



d 3 q 



1 f 9 ^ d 6 i 

^ (P) = E p (4E$ -Ml) I q^^ (P + q) W 



(A13) 



In the non relativistic limit one has 



M 2 D = (2m - B) 2 ~ 4m(m - B), E p = (p 2 + m 2 ) 1 / 2 ~ m 



(A14) 



and one obtains 



(^ +fl ) w( p) = ^y 



<? 2 / . \ A 
2<^p(p + q) 



q 2 +{!■ 



(2vr) 3 



(A15) 



which may be written as 



p z f d 3 q 

— ¥>p(p) + J V r (q)(^ P (p + q)^3 = E(p P (p) 



2m R 



(A16) 



where = m/2 is the reduced mass of the two nucleons, E = —B and V(q) is obtained by 
comparing the two equations above. (fp(p) is hence the bound state solution of the Schrodinger 
equation of the two nucleons' system. One may notice that, as in the scalar case considered here 
<!>£> has the dimensions of a mass, $p has dimensions mass -3 and ipp mass -2 . The bound state 
solution of the Schrodinger equation of the two nucleons' system <fp(p) cannot hence be identified 
with the usual non relativistic deuteron wave function normalized to one after integration over p 
(which has dimensions mass -3 / 2 ). Obviously in the non relativistic limit, for any a, any function 
Epipp(p) becomes a solution of the Schrodinger equation of the two nucleons' system. The 
connection with the usual non relativistic deuteron wave function may hence be obtained by 
determining the value of a through the normalization condition of the deuteron wave function: 
The link between ip p and the vertex function is given by 



*d(p) = 



— + p ) — m + it 



P 

~2 



P 



ig fI3 . , d A q 

— 2 X{P, P + Q) 17, Y4 

q A — n z (27r) 



= -iE p {AE 2 -M 2 D )<p P (p) 
in such a way that one may write 



x(p,p) 



— m + it 

-id 2 , \ 
q 2 + p z (2tt) 6 



d?q 



(A17) 



21 



<Mp) 



**r>(p) 



(A18) 



£p(4£|-M 2 ) 

keeping into account the normalization of the Bethe-Salpeter wave function, mamely the relation 



**D 



(M 2 D -4E 2 ) 2 \2E p (2^ 



d 3 p 



= 1 



(A19) 



one can make the position 



<Mp 2 ) 



E p {AEj-Ml) (2tt)3/2 



(A20) 



where <fp(p 2 ) is hence normalized to 1 and, in the non relativistic limit, it is a solution of the 
Schrodinger equation of the two nucleons' system. The solution of the Bethe-Salpeter equation 
may hence be expressed in terms of <fp{p 2 ) as 



*p(p) 



(?-')■ 



to 2 + ie 



(A21) 



while the function ^p,(p), defined in Eq.5, 6 is expressed in terms of ipp(p 2 ) by the relation 



jMg) 

p- 



j [D-pf- 



(15 a 



1 



m 2 + ie] [p 2 — m 2 + ie] 2ir p_ (D — p) 2 — m 2 



(A22) 



where p is the three momentum in the deuteron rest frame and E p = \Jp 2 + m 2 . 

The three momentum in the deuteron center of mass frame can be expressed covariantly. In 
the case of interest one of the two nucleons is on shell. If, as shown in Fig. 8, one puts t = {D—p) 2 
one obtains 



2 , 2 

p + m = 



4M 2 



{t-m 2 - M 2 D 



2 \2 



(A23) 



Introducing light cone variables, with Z the momentum fraction of the virtual nucleon with 
respect to half of the deuteron momentum, one obtains 



—^-(—Mf) + — m? 
4M 2 V2 D Z K 



m 



(A24) 



and, for E p and p z 
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virtual 



real 



Figure 8: Deuteron vertex, one nucleon is virtual and one nucleon is on shell 



, , {—Mi mj 

2M D \2 u Z l ) 1 2M D \2 D Z 1 



Pz 



(A25) 



The function ^ d (p) is hence finally expressed in terms of light cone variables through the non 
relativistic deuteron wave function <^p(p 2 ) as 



m D {zD_- Pt ) = (2vr) 



3/2 



1 /Z 



2M D 



x cp P 



(y4M 2 D ( 



z mt 

2 



m 



(A26) 



Notice that, since in principle ^£>(Z, B) = 'I'd (2 — Z, B), the expression above has to be under- 
stood as symmetrized for the exchange Z 44> (2 — Z). 
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